METHOD AND UNIT 
FOR 

SUPPRESSING A PERIODIC PATTERN 



BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a periodic - pattern 
suppression method and unit, and more particularly to a method and 
unit for reducing and removing a stripe pattern (including a moire 
component) , corresponding to a stationary grid, from an image 
photographed, for example, by use of the stationary grid, the image 
having the stripe pattern superposed on the image of a subject. 

Description of the Related Art 

Radiation recording - reproducing systems have been 
proposed in Japanese Unexamined Patent Publication Nos . 
55 (1980) -12429, 56 (19 81) -11395, 55 (1980) -163472, 56 (1981)- 
164645, 55 (1980) - 11634 0 , etc. These systems utilize a 
storage- type phosphor ( s timulatable phosphor) , which stores part 
of radiation energy when irradiated with radiation, such as X- 
rays , etc. , and emits photos timula ted luminescent light according 
to the amount of the stored radiation energy when irradiated with 
excitation light such as visible light, etc. The radiation image 
of a subject, such as a human body, etc., is temporarily 
photographed and recorded on a storage- type fluorescent sheet. 
Then, excitation light, such as laser light, etc., is scanned on 
the storage-type fluorescent sheet to obtain photos timulated 



luminescent light. The obtained photos timulated luminescent 
light is photoelec tr ically read out and converted to an electrical 
image signal by a reading means such as a photomul tipl ier , etc. 
Based on the image signal, the radiation image of the subject is 
output as a visible image to a recording material such as a 
photosensitive material, etc., or to a cathode-ray tube (CRT) 
display unit, etc. 

There are cases where, when photographing and recording 
the radiation image of a subject on the above-mentioned 

hi- :i 

^^^lO storage- type fluorescent sheet, photographing is performed with 

l"IJ 

a stationary grid disposed between the subject and the fluorescent 
^ sheet. In the stationary grid, lead plates, which do not transmit 

jj^ radiation, and aluminum or wood plates, which easily transmit 

H radiation, are alternately disposed at fine pitches of about 4 

15 plate/mm so that radiation scattered by the subject is not 

irradiated to the fluorescent sheet. If photographing is 
performed using this stationary grid, radiation scattered by the 
subject will become less liable to be irradiated to the fluorescent 
sheet and therefore the contrast of the radiation image of the 
20 subject can be enhanced. However, the grid image in the form of 

a fine stripe pattern corresponding to the stationary grid, along 
with the subject image, is recorded. 

Because of this, the applicant of this application has 
proposed, in Japanese Unexamined Patent Publication No. 
25 10 (1998) - 164737 , a method of obtaining a s tripe - reduced image 

which is easy to observe, by applying a filtering process for 
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removing a spatial frequency component corresponding to the stripe 
pattern of a stationary grid. In this method, for example, in the 
case where the grid array pitch (distance between the centers of 
adjacent grid rods) of a stationary grid (hereinafter referred to 
as a grid pitch) is 4 grid - 1 ine/mm, a stripe pattern develops in 
a spatial frequency band near 4.0 cycle/mm. In order to remove 
the stripe pattern, a filtering process is performed by a filter 
which removes or reduces its response at this frequency band. 

However, the filter employed in the method in the 



t-ff 

Quo aforementioned publication No. 10 (1998) - 164737 suppresses and 

removes not only the stripe pattern of the stationary grid but also 
the high frequency component contained in the original image, 
because its response at a frequency component near a spatial 
frequency corresponding to the grid pitch of a stationary grid used, 
15 and its response at a high frequency component greater than that, 

are made zero. Thus, there is a problem that an image will be 
reduced in sharpness. 

SUMMARY OF THE INVENTION 
The present invention has been made in view of the 
20 aforementioned circumstances. Accordingly, it is the primary 

object to provide a periodic - pa ttern suppression method and a 
periodic -pattern suppression unit which are capable of making a 
periodic pattern (such as a stripe pattern, which develops in an 
image because of a stationary grid, etc., and the like) 
25 inconspicuous without reducing sharpness. 

To achieve this end and in accordance with a first 
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important aspect of the present invention, there is provided a 
periodic -pattern suppression method of reducing a spatial 
frequency component which forms a periodic pattern contained in 
an original image signal, the method comprising the steps of: 

transforming the original image signal, represented in 
a real space domain, into a plurality of transformed image signals 
which can be handled in a frequency domain; and 

reducing a transformed image signal of the transformed 

5., 3. 

11 image signals which has a desired frequency range containing a 

lUO spatial frequency component corresponding to at least a frequency 

Ij 

^ of the periodic pattern in only the vicinity of an array direction 

J of the periodic pattern. 

y The words "periodic pattern" are intended to broadly mean 

i t 

3 patterns which develop periodically in an image. For example, the 

15 periodic pattern indicates a great variety of patterns such as a 

stripe pattern, which occurs in an image because of a stationary 
grid, a moire image, which occurs in an image because of the 
characteristics of an imaging system in a television set, and the 
like . 

20 In accordance with a second important aspect of the 

present invention, there is provided a periodic -pattern 
suppression method of reducing a spatial frequency component 
resulting from a stationary grid, contained in an original image 
signal photographed using the stationary grid, the method 

25 comprising the steps of: 

transforming the original image signal, represented in 
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a real space domain, into a plurality of transformed image signals 
which can be handled in a frequency domain; and 

reducing a transformed image signal of the transformed 
image signals which has a desired frequency range containing a 
spatial frequency component corresponding to at least a grid array 
frequency of the stationary grid in only the vicinity of a grid 
array direction of the stationary grid. 

The expression "transforming into a plurality of 
transformed image signals" means transforming into signals to 
0*110 which a desired process for a desired frequency component contained 

I y 

H= in the original image signal can be applied. Particularly, in the 

i present invention, it means transforming into signals on which a 

process of suppressing a frequency component, contained in the 
original image signal represented in a real space domain, which 
15 corresponds to the stationary grid can be applied. For example, 

an original image signal expressed in a real space domain can be 
transformed into image signals expressed in a frequency domain 
(frequency spectra) by Fourier transformation. Also, an original 
image signal expressed in a real space domain can be transformed 
20 into subband signals expressed in a frequency domain by 

mul tiresolution transformation (employing a wavelet transform or 
a Laplacian pyramid expansion) . 

The expression "spatial frequency component resulting 
from a stationary grid" includes not only the spatial frequency 
25 component of the stationary grid itself, but also includes a moire 

component which occurs due to the stationary grid by sampling at 
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sampling cycles less than a Nyquist frequency or by a reducing 
process. This is why the expression "spatial frequency component 
corresponding to a grid array frequency" is employed. This 
expression means is intended to include not only the same spatial 
5 frequency component as the grid frequency but also the same spatial 

frequency component as the moire frequency related to this. 
_ The expression "reducing a transformed image signal of 

y 

the transformed image signals which has a desired frequency range 
^2 containing a spatial frequency component corresponding to at least 

5l0 a grid array frequency of the stationary grid" means to apply a 

flJ 

H" process for suppressing both a spatial frequency component 

B corresponding to the grid frequency and a spatial frequency 

ry component near the grid frequency, and does not mean to suppress 

m 

y almost all high frequency components greater than a spatial 

15 frequency component corresponding to the grid frequency, as in the 

conventional method. That is, in the present invention, as fas 
as it is possible, high frequency components (e.g., a spatial 
frequency component corresponding to a Nyquist frequency and 
spatial frequency components near the Nyquist frequency) greater 
20 than a spatial frequency component corresponding to the grid 

frequency are not suppressed. 

The expression "in only the vicinity of a grid array 
direction of the stationary grid" means to apply the suppressing 
process in the grid array direction of the stationary grid or the 
25 neighboring directions, and means not to apply the suppressing 

process in the directions other than that, unlike applying the 
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suppressing process independently of the grid direction, as in the 
conventional method. The "stationary grid" used herein, as 
described later, may be a stationary grid which is actually used, 
or a stationary grid, differing in grid direction, which is to be 
5 used. Therefore, these stationary grids are included in the 

"stationary grid" in the expression "in only the vicinity of a grid 
array direction of the stationary grid." In other words, since 
^13 the grid direction of a stationary grid that is used in ordinary 

gn units has been determined to some degree and is, for example, a 

in 

yfUO horizontal or vertical scanning direction, the process of 

M suppressing a component corresponding to the grid frequency may 

□ be applied in the horizontal scanning direction and the vertical 

fy scanning direction orthogonal to each other. To speak in plainer 

fU 

O language, the suppressing process may be applied in almost all grid 

15 directions of a stationary grid to be used. For instance, in 

Fourier space with a horizontal scanning direction as a v-axis and 
a vertical scanning direction as a u-axis, only a desired spatial 
frequency component containing a grid frequency in the vicinity 
of the v-axis (having a slight width in the positive and negative 
20 directions of the u-axis orthogonal to the v-axis) and in the 

direction of the v-axis, and a desired spatial frequency component 
containing a grid frequency in the vicinity of the u-axis (having 
a slight width in the positive and negative directions of the v-axis 
orthogonal to the u-axis) and in the direction of the u-axis, may 
25 be suppressed. In this case, the grid length direction of one 

stationary grid is the same as the grid array direction of the other 
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stationary grid, and consequently, the suppressing process is 
applied in both of the grid array direction and grid length 
direction of a stationary grid orthogonal to each other. 

In the periodic - pattern suppression method according to 
the present invention, it is preferable to perform the suppressing 
process only on a stationary grid that is actually used. 
Performing the suppressing process only on a stationary grid that 
is actually used shall mean to apply the suppressing process only 
in the grid array direction of a stationary grid, which is actually 
used, or the neighboring directions, and shall mean not to apply 
the suppressing process in the directions (including also the grid 
length direction of the stationary grid) other than that. For 
instance, in the case of a vertical stationary grid, it means to 
suppress only a component which has a desired spatial frequency 
range containing a grid frequency in the vicinity of the v-axis 
and in the direction of the v-axis. 

In a preferred form of the periodic -pattern suppression 
method according to the present invention, the aforementioned 
transforming step obtains a plurality of transformed image signals 
by applying two-dimensional wavelet transformation to the original 
image signal by the use of a low-pass filter which splits a band 
so that its response at a frequency greater than the spatial 
frequency of the stationary grid becomes approximately zero . Also, 
the aforementioned reducing step further applies a process of 
reducing a component less than a predetermined frequency and then 
performs inverse wavelet transformation, with respect to a signal 



8 




of the transformed image signals which contains a spatial frequency 
component corresponding to the grid array frequency. 

In this case, it is more desirable that the 
aforementioned reducing step reduce a component less than the 
5 predetermined frequency, by recursively and repeatedly applying 

one - dimensional wavelet transformation to the transformed image 
signal (wavelet transform coefficient signal) , containing a 
v3 spatial frequency component corresponding to the grid array 

ui frequency, in the grid length direction of the stationary grid by 

pIlO a predetermined number of times by the use of a predetermined band 

M= splitting filter, then making zero transform coefficients of a low 

O frequency image signal of a plurality of image signals obtained 

flj by one- dimensional wavelet transformation, and applying inverse 

Q one - dimensional wavelet transformation. 

5 : 

15 The "predetermined band splitting filter" does not 

always need to be a low-pass filter which splits a band so that 
its response at a frequency greater than the spatial frequency of 
a stationary grid becomes approximately zero. 

In applying one - dimensional wavelet transformation in 

20 the grid length direction of the stationary grid, when the 

stationary grid is a vertical grid a HL component obtained by 
two-dimensional wavelet transformation is subjected to one - 
dimensional wavelet transformation in the vertical scanning 
direction, and when it is a horizontal grid an LH component obtained 

25 by two-dimensional wavelet transformation is subjected to one- 

dimensional wavelet transformation in the horizontal scanning 
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direction. On the other hand, in the case of a cross grid, a HH 
component obtained by two-dimensional wavelet transformation is 
subjected to one - dimensional wavelet transformation in the 
horizontal scanning direction and is then subjected to one- 
5 dimensional wavelet transformation in the vertical scanning 

direction. Note that the order of the horizontal and vertical 
scanning directions may be reversed. 

u 

iQ Besides such a method, a component less than a 

pi predetermined frequency representing the stationary grid 

OHO component can also be reduced by extracting only a component 

nj 

representing an original image from an image signal, which contains 
p the stationary grid component, among the image signals obtained 

J. — 

ilJ by wavelet transformation. 

RJ 

p In a preferred form of the periodic -pattern suppression 

i — ■- 

15 method according to the present invention, the reducing step 

calculates powers of a plurality of transformed image signals, 
judges the grid length direction of the stationary grid, based on 
whether or not each of the calculated powers is greater than a 
predetermined threshold value, and applies the process for 

20 reducing a component less than a predetermined frequency, based 

on the result of judgement. 

The periodic -pattern suppression method according to 
the present invention is not limited to a stationary grid that is 
actually used. For each grid which is to be used, the 

25 aforementioned reducing step may be performed by employing the 

transformed image signals obtained by the aforementioned two- 
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dimensional wavelet transformation . 

In the periodic -pa ttern suppression method according to 
the present invention, the aforementioned transforming step may- 
obtain a plurality of transformed image signals by applying 
5 one - dimensional wavelet transformation to the original image 

signal in the grid length direction of the stationary grid by the 
use of a predetermined band splitting filter. Also, the 
i5 aforementioned reducing step may further apply a process of 

gi reducing a component less than a predetermined frequency and then 

jlllO perform inverse wavelet transformation, with respect to a low 

3 y 

tj, frequency image signal of the transformed image signals which 

a 

rj contains a spatial frequency component corresponding to the grid 

h[j array frequency of the stationary grid. 

nil ..... 

similarly, the above case is not limited to a stationary 

hi 

15 grid that is actually used. The aforementioned transforming step 

and the reducing step can be performed on each stationary grid to 
be used. 

In accordance with a third important aspect of the 
present invention, there is provided a per iodic - pat tern 
20 suppression unit for reducing a spatial frequency component which 

forms a periodic pattern contained in an original image signal, 
the unit comprising the steps of: 

image signal transformation means for transforming the 
'original image signal, represented in a real space domain, into 
25 a plurality of transformed image signals which can be handled in 

a frequency domain; and 
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per iodic - pattern - component suppression means for 
reducing a transformed image signal of the transformed image 
signals which has a desired frequency range containing a spatial 
frequency component corresponding to at least frequency of the 
5 periodic pattern in only the vicinity of an array direction of the 

periodic pattern. 

In accordance with a fourth important aspect of the 
5 present invention, there is provided a per iodic -pattern 

m suppression unit for reducing a spatial frequency component 

CfilO resulting from a stationary grid, contained in an original image 

:H 5 

tl signal photographed using the stationary grid, the unit 

pj comprising: 

kj image signal transforming means for transforming the 

flj 

^ original image signal, represented in a real space domain, into 

15 a plurality of transformed image signals which can be handled in 

a frequency domain; and 

stationary - grid suppressing means for reducing a 
transformed image signal of the transformed image signals which 
has a desired frequency range containing a spatial frequency 
20 component corresponding to at least a grid array frequency of the 

stationary grid in only the vicinity of a grid array direction of 
the stationary grid. 

In a preferred form of the periodic -pattern suppression 
unit according to the present invention, the image signal 
25 transforming means obtains the plurality of transformed image 

signals by applying two-dimensional wavelet transformation to the 
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original image signal by the use of a low-pass filter which splits 
a band so that its response at a frequency greater than the spatial 
frequency of the stationary grid becomes approximately zero . Also , 
the stationary - grid suppressing means further applies a process 
5 for reducing a component less than a predetermined frequency and 

then performs inverse wavelet transformation, with respect to an 
image signal of the transformed image signals which contains a 

LJ 

43 spatial frequency component corresponding to the grid array 

g1 frequency of the stationary grid. 

W . . . 

ijlio In this case, it is desirable that the s tat lonary - grid 

ni 

suppressing means reduce a component less than the predetermined 
O frequency, by recursively and repeatedly applying one - dimensional 

nj wavelet transformation to the transformed image signal, containing 

a spatial frequency component corresponding to the grid array 
15 frequency, in a grid array direction of the stationary grid by a 

predetermined number of times by the use of a predetermined band 
splitting filter, then making zero transform coefficients of a low 
frequency image signal of a plurality of image signals obtained 
by the one - dimensional wavelet transformation, and applying 
20 inverse one - dimensional wavelet transformation. 

The periodic - pattern suppression unit according to the 
present invention may further include stationary grid - direction 
judging means for calculating powers of the plurality of 
transformed image signals and judging the grid length direction 
25 of the stationary grid, based on whether or not each of the 

calculated powers is greater than a predetermined threshold value. 
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In this case, the stationary grid - direction judging means applies 
the process for reducing a component less than a predetermined 
frequency, based on the judgement made by the stationary grid- 
direction judging means. 
5 The above-mentioned stationary - grid suppressing means, 

in addition to a stationary grid which is actually used, is capable 
of applying the process for reducing a component less than a 
yi predetermined frequency, to each stationary grid to be used. 

•IX 

Oil In still another preferred form of the periodic - pa t tern 

' ''\ s 

•7:1' i 

iJilO suppression unit according to the present invention, the image 

M signal transforming means obtains the plurality of transformed 

O image signals by applying one - dimensional wavelet transformation 

fiy to the original image signal in the grid length direction of the 

p stationary grid by the use of a predetermined band splitting filter . 

15 Also, the stationary -grid suppressing means further applies a 

process for reducing a component less than a predetermined 
frequency and then performs inverse wavelet transformation, with 
respect to a low frequency image signal of the transformed image 
signals which contains a spatial frequency component corresponding 
20 to the grid array frequency of the stationary grid. 

In the above case, the image signal transforming 
means may apply the one - dimensional wavelet transformation in the 
grid length direction of each stationary grid to be used, and the 
stationary-grid suppressing means may apply the reducing process 
25 and the inverse wavelet transformation to each stationary grid to 

be used. 
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According to the periodic -pattern suppression method 
and unit of the present invention, an image signal which represents 
an image having a periodic pattern is converted into a plurality 
of transformed image signals which can be handled in a frequency 
5 domain. Also, a transformed image signal of the transformed image 

signals, which has a desired frequency range containing a spatial 
frequency component corresponding to at least a frequency of the 
a5 periodic pattern in only the vicinity of an array direction of the 

pi periodic pattern, is reduced. As a result, the periodic pattern 

m 

gilO in an image can be made inconspicuous. 

In addition, in the case where a periodic pattern results 

Q from a stationary grid, an original image signal, photographed 

p™ 

f|3 using a stationary grid, and represented in a real space domain, 

fy 

is transformed into a plurality of image signals which can be 
15 handled in a frequency domain. Also, among the plurality of image 

signals, the image signal, which has a desired frequency range 
containing the same spatial frequency component as the grid 
frequency of the stationary grid in approximately the grid array 
direction of the stationary grid, is reduced. Therefore, a stripe 
20 pattern resulting from the stationary grid will become 

inconspicuous. On the other hand, because all frequency 
components greater than the spatial frequency component resulting 
from the stationary grid are not suppressed, by that amount an image 
signal finally obtained will contain components near the spatial 
25 frequency component corresponding to a Nyquist frequency and 

therefore there is no possibility that sharpness will be reduced. 
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Unlike performing the suppressing process in all 
directions on Fourier space independently of grid direction, as 
in the conventional method, the stationary grid component is 
removed only in the grid array direction of the stationary grid 
and therefore there is no possibility that the image will be 
deteriorated more than necessary. 

In the case where grid direction is known, there is no 
possibility that that the image will be deteriorated even more than 
necessary, if the aforementioned suppressing process is given only 
in the grid array direction of the stationary grid. 

If the suppressing process is applied in the grid array 
directions of stationary grids (which differ in grid direction) 
to be used, the effect of the aforementioned suppressing process 
can be obtained without taking into consideration the grid 
direction of a stationary grid which is actually used. 
Particularly, the effect is great when two-dimensional wavelet 
transformation is performed. 

If an image signal is given two-dimensional wavelet 
transformation by the use of a low-pass filter which splits a band 
so that its response at a frequency greater than the spatial 
frequency of the stationary grid becomes approximately zero, and 
if, among the image signals obtained by the wavelet transformation, 
the image signal containing a stationary grid component is further 
subjected to a process of reducing a component less than a 
predetermined frequency, an image with a reduced stationary grid 
component can be restored when the original image (resolution level 




0) is restored by the use of the signal in which the component less 
than the predetermined frequency has been reduced. 

Since the low - resolution signal component obtained by 
the wavelet transformation at the initial stage does not contain 
5 the stationary frequency component of the stationary grid, each 

wavelet - transform coefficient signal will not contain the spatial 
frequency component of the stationary grid, even if low - resolution 
y;5 image signals are obtained by further applying wavelet 

Bl transformation over a plurality of stages. Therefore, even if an 

Olio image is restored to an intermediate resolution level without being 

flj 

M: restored to resolution level 0 when it is restored, there is no 

p possibility that a moire pattern resulting from the stationary grid 

fU will occur in a reduced - scale image ( low- resolution image), 

p because no grid component is contained in low - resolution signals 

15 of level 1 and levels thereafter. That is, even if an image is 

arbitrarily enlarged or reduced, no moire pattern will occur. 

If, among a plurality of image signals obtained by 
recursively and repeatedly applying one - dimensional wavelet 
transformation to the image signal (wavelet - transform coefficient 
20 signal), containing the stationary grid component, in the grid 

direction by a predetermined number of times by the use of a 
predetermined band splitting filter, the transform coefficients 
of the low frequency image signal are made zero, in the case where 
image signals that can be handled in a frequency domain are obtained 
25 by wavelet transformation, a restored image can be made an image 

with only the stationary grid component suppressed, as the number 
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of repeats is increased. Thus, an image with even higher sharpness 
can be obtained. 

If powers of a plurality of image signals obtained by 
wavelet transformation are calculated and the grid direction is 
judged based on whether or not each of the calculated powers is 
greater than a predetermined threshold value, when the suppressing 
process is performed it is not necessary to previously know what 
kind of stationary grid is used in photographing. 

Furthermore, if an image signal, photographed using a 
WHO stationary grid, and represented in a real space domain, is 

subjected to one - dimensional transformation in the grid length 
direction of the stationary grid by the use of a predetermined band 
splitting filter, and if, among a plurality of image signals 
obtained by the wavelet transformation, the image signal on the 
15 low frequency side containing the stationary grid component is 

further subjected to the process of reducing a component less than 
a predetermined frequency, an image signal which has a desired 
frequency range containing a spatial frequency component resulting 
from the stationary grid can be reduced. When an image is restored 
2 0 by the use of the signal in which the component resulting from this 

stationary grid has been reduced, an image with a reduced grid 
component can be restored regardless of resolution level at which 
the image is restored. 

BRIEF DESCRIPTION OF THE DRAWINGS 
25 The present invention will be described in further detail 

with reference to the accompanying drawings wherein: 



18 



FIG. 1 is a top view showing a radiation image 
photographing apparatus ; 

FIG. 2 is a diagram showing the image of a subject and 
the image of a grid photographed by the radiation image 
photographing apparatus shown in FIG. 1; 

FIG. 3 is a perspective view showing a radiation image 
reader to which a periodic -pa ttern suppression unit according to 
the present invention, for carrying out a per iodic - pa ttern 
suppression method, is applied; 

FIG. 4 is a diagram showing the relationship between 
scanning directions and an image to be read; 

FIG. 5 is a block diagram showing an example of the image 
signal processing section of the per iodic - pa ttern suppression unit 
of the present invention; 

FIG. 6 is a block diagram showing how a wavelet transform 
process is performed by a wavelet transform section; 

FIG. 7 is a block diagram showing wavelet transform 

means ; 

FIG. 8 is a graph showing the frequency response 
characteristics of the low-pass filters in wavelet transform 
filters ; 

FIG. 9A is a diagram showing the decomposed components 
of the original image signal obtained after the original image 
signal has been subjected to the first two-dimensional wavelet 
transformation ; 

FIG. 9B is a diagram showing the state in which the 



original image signal has further been decomposed from the state 
of FIG. 9A; 

FIG. 10 is a schematic diagram showing how the effect 
of one - dimensional wavelet transformation is obtained by repeating 
5 the transformation; 

FIG. 11 is a block diagram showing the construction of 
an inverse wavelet transform section; 
yp FIG. 12 is a block diagram showing how inverse wavelet 

gi transformation is performed by the wavelet transform section shown 

ij^lO in FIG. 11; 

fU 

i^. FIG. 13 is a diagram showing advantages of the present 

ri invention; and 

m FIG. 14 is a diagram showing another embodiment of the 

nj 

=-1 present invention. 

U 

15 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A preferred embodiment of the present invention will 
hereinafter be described in detail with reference to the drawings. 
The following description will be made according to a mode, in which 
in a radiation image information recording - reproducing system 

20 utilizing a storage-type fluorescent sheet as a recording sheet, 

the radiation image of a human body recorded on the recording sheet 
is read out as a digital image signal by laser-beam scanning, as 
described in the aforementioned Japanese Unexamined Patent 
Publication Nos . 55 ( 1980) - 12429 , 56 ( 19 8 1 ) - 1 13 9 5 , etc. 

25 Fig. 1 schematically shows a radiation image 

photographing apparatus. As shown in the figure, radiation 2 
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emitted from a radiation source 1 passes through a subject 3 and 
reaches a grid 4. In the grid 4, lead plates 4a which absorb the 
radiation 2 and aluminum plates 4b which transmit the radiation 
2 are alternately disposed at pitches of 4 plate/mm with a slight 
tilt in accordance with position so that the radiation 2 emitted 
from the radiation source 1 is incident straight on a recording 
sheet 11 via the aluminum plates 4b (refer to Fig. 1) . Because 
of this, the radiation 2, emitted from the radiation source 1, and 
Oil transmitted straight through the subj ect 3 , is absorbed and blocked 

Olio by the lead plates 4, but it is transmitted through the aluminum 

M= plates 4b and irradiated to the recording sheet 11. As a result, 

O a grid image 6 with four stripes per millimeter, along with the 

image of the subject 3, is stored and recorded on the recording 
sheet 11. On the other hand, radiation 2a scattered within the 
15 subject 3 is incident obliquely with respect to the tilt of the 

grid 4. Since the scattered radiation 2a incident on the aluminum 
plates 4b is absorbed in the lead plates 4a within the grid 4, or 
it is reflected at the surface of the grid 4, it is not irradiated 
to the recording sheet 4. Thus, a distinct radiation image with 
20 less irradiation of the scattered radiation 2a is stored and 

recorded on the recording sheet 11, 

Fig. 2 illustrates a subject image (indicated by oblique 
lines) 5 and a grid image 6 in the form of a stripe pattern, stored 
and recorded on the recording sheet 11 by performing photographing, 
25 using the grid 4. Thus, a radiation image with the subject image 

5 superposed on the grid image 6 is recorded on the recording sheet 
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11 . 

Fig. 3 illustrates a radiation image reader to which an 
embodiment of a periodic - pat tern suppression unit according to the 
present invention, for carrying out a periodic -pattern suppression 
5 method, is applied. 

The recording sheet 11 with the radiation image, set at 
a predetermined position in a reading section 10, is conveyed in 
iO a direction of arrow Y at intervals of 10 scan/mm by sheet conveying 

01 mean 10, such as an endless belt, which is driven by drive means 

yllO (not shown) . On the other hand, a light beam 17 emitted from a 

ki. laser light source 16 is reflected and deflected by a rotating 

Q polygon mirror 18 which is driven to rotate at high speeds in an 

ry arrow - indicating direction by a motor 24 . The deflected light beam 

IIJ 

Q 17 is passed through a focusing lens 19 such as fO , for example. 

15 The light beam 17 is changed in direction by a mirror 20 and is 

incident on the sheet 11 and moved in a direction approximately 
perpendicular to the vertical scanning direction (direction of Y 
arrow) . From a position on the sheet 11 irradiated with the light 
beam 17, photos timulated luminescent light 21 with a quantity of 

20 light corresponding to the radiation image information being 

stored and recorded is emitted and incident on the incidence end 
22a of an optical guide 22. The photos timulated luminescent light 

21 is totally reflected at the internal surface of the optical guide 

22 and arrives at the emergence end 22b of the optical guide 22 
25 and is received by a photomul t ipl ier 23. The photos timulated 

luminescent light 21 representing the radiation image is 
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photoelectrically detected and converted to an electrical signal 

The output analog signal is logarithmically amplified 
by a log amplifier 26. The amplified signal is sampled and 
digitized at sampling intervals corresponding to a spatial 
frequency of fs = 10.0 cycle/mm by an A/D converter 28. In this 
manner, a digital image signal S carrying a high-density image with 
a high reading density is obtained. Note that the image signal 
S contains the information of a spatial frequency band of 4.0 
cycle/mm (which is information on the grid image 6 shown in Fig. 
2) lower than the highest spatial frequency (Nyquist frequency to 
be described later) f n = 5.0 cycle/mm in a desired spatial frequency 
range, required to reproduce and output a satisfactory visible 
radiation image. Also, the information on the grid image 6 should 
be removed because it is one of the causes making it difficult to 
view the visible radiation image when observing it. 

The image signal S represents image information obtained 
by scanning the laser beam on the sheet 11 two-dimensionally , as 
shown in Fig. 4. That is, the sheet 11 is moved in the vertical 
scanning direction (direction of arrow Y) , while the laser beam 
is being scanned on the sheet 11 in the horizontal scanning 
direction (direction of arrow X) . The image signal S thus obtained 
carries information less than Nyquist frequency fn, so it also 
contains the information (4.0 cycle/mm) on the grid image 6 shown 
in Fig. 2. In this embodiment, the moire of the grid image 6 due 
to aliasing will not occur because the analog signal is digitized 
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at sampling intervals corresponding to the spatial frequency fs 
which is two or more times the frequency (4.0 cycle/mm) of the 
information on the grid image 6 . 

After being stored temporarily in a storage section 29, 
the image signal S is input to an image signal processing section 
30 and processed in the following manner. 

Fig. 5 shows the image signal processing section (one 
form of a periodic -pattern suppression unit) 30 for carrying out 
the per iodic -pattern suppression. method of the present invention. 
filO As shown in the figure, the image signal processing section 30 has 

(1) mul tiresolut ion decomposing means 32 as image signal 
□ transforming means for applying a mul tiresolut ion decomposing 

process to the image signal S read out from the storage section 
29, (2) processing means 33 for applying a desired process (e.g., 
15 an enhancing process) to an image signal which has a predetermined 

frequency band (band - 1 imi ted image signal) , among the band- limited 
image signals (band-pass and sub-band signals) decomposed to a 
plurality of frequency bands by the mul tiresolution decomposing 
means 32, (3) restoring means 34 for obtaining a processed image 
20 signal by restoring both the band-limited image signal, subjected 

to the desired process by the processing means , and the band - limited 
image signals carrying other frequency bands , and (4) image output 
means 35 for reproducing a visible image, based on the processed 
image signal restored by the restoring means 34. 
25 Note that in this embodiment, the mul tiresolution 

decomposing process employing a wavelet transform process is used 
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when the image signal S represented in an actual spatial domain 
is transformed into image signals which can be handled in a 
frequency domain. Therefore, the mul tiresolution decomposing 
means 32 functions as wavelet transform means for applying wavelet 
transformation to the image signal S representing the radiation 
image, while the restoring means 34 functions inverse wavelet 
transform means for obtaining a processed image signal by applying 
an inverse wavelet transform process. In the following 
description, the mul tiresolution decomposing means 32 is referred 
to as a wavelet transform section 32 and the restoring means 34 
is referred to as an inverse wavelet transform section 34 . 

Fig. 6 shows how the wavelet transform process as the 
mul tiresolution decomposing process is performed by the wavelet 
transform section 32. In Fig. 7, each wavelet transform means 32a 
is shown in detail. As shown in Fig. 6, the wavelet transform 
section 32 is provided with a plurality of wavelet transform means 
32a corresponding in number to resolution (pixel density) levels. 
Similarly, the inverse wavelet transform section 34 is provided 
with a plurality of inverse wavelet transform means 34a 
corresponding in number to the plurality of wavelet transform means 
32a . 

As also shown in Fig. 6, the image signal processing 
section 30 is equipped with stationary grid - direction judging 
means (hereinafter referred to as direction judging means) 35 and 
stationary grid - component suppressing means 36. The direction 
judging means 35 is used for calculating powers of signals HLl, 



LHl, and HHl obtained by wavelet transformation, and for judging 
the grid length direction of the stationary grid 4 ( i.e., grid 
direction) , based on whether each of the calculated powers is 
greater than a predetermined threshold value HTl . The stationary 
grid - component suppressing means 36 is used to reduce an image 
signal of the image signals HL, LH , and HH which has a desired 
frequency range containing at least a spatial frequency component 
resulting from the stationary grid 4 in the grid array direction 
of the stationary grid 4 (i.e., a direction perpendicular to the 
grid direction) . 

Now, the periodic -pattern suppression method according 
to the present invention will be described in detail. This 
embodiment performs two-dimensional wavelet transformation in 
which wavelet transform coefficients are orthogonal to one 
another . 

As shown in Figs. 6 and 7, if the digital image signal 
S, read out from the storage section 29, which represents the 
original image, is input as an original image signal 3^^,^ to the 
wavelet transform section 32, the original image signal S^^^ is 
subjected to wavelet transformation. More specifically, a 
filtering process is performed in the horizontal scanning 
direction of the original image signal S^^i (equivalent to a signal 
LLO) by wavelet functions HI, Gl, and the pixels in the horizontal 
scanning direction are thinned out at every other pixel, whereby 
the number of pixels in the horizontal scanning direction is reduced 
to one half. Note that the thinned-out signals obtained in this 



way are represented by "t2" in Fig. 7. The function HI is a 
high-pass filter, while the function Gl is a low-pass filter. 
Furthermore, the filtering process is performed on each of the 
thinned-out signals in the vertical scanning direction by the 
5 wavelet functions HI, Gl , and the pixels in the vertical scanning 

direction are thinned out at every other pixel, whereby the number 
of pixels in the vertical scanning direction is reduced to one half . 
y3 In this manner, wavelet - transform coefficient signals 

01 (hereinafter also referred to simply as signals) HHl, HLl, LHl, 

ffUO and LLl are obtained. The signal LLl represents a one -fourth 

y= reduction image obtained by reducing the original image to one half 

O in both the vertical and horizontal directions . In the respective 

|1J one- fourth reduction images of the original image, the signal LHl 

Q represents an image of a high frequency component (horizontal edge) 

15 in the vertical scanning direction, the signal HLl represents an 

image of a high frequency component (vertical edge) in the 
horizontal scanning direction, the signal HHl represents an image 
of a high frequency component (diagonal edge) in the diagonal 
direction, and the signal LLl represents an image of a low frequency 
20 component of half resolution with respect to the original image. 

The frequency which becomes a reference for band split, i.e., the 
frequency at the boundary between the low frequency component, and 
each of the high frequency components is determined by the filter 
characteristics of the functions HI, Gl . In the filtering process 
25 by the functions HI, Gl at the initial stage, the function Gl as 

a low-pass filter has such a characteristic that its response at 
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a spatial frequency of 4.0 cycle/mm or more is practically zero 
at the grid pitches of the stationary grid 4. Also, the function 
HI as a high-pass filter has a high-pass characteristic for 
compensating for the low-pass characteristic of the function Gl . 
5 As a filter with such a characteristic, there is, for example, a 

filter having the same characteristic as a moire removing filter 
described in Japanese Unexamined Patent Publication No. 

f'l 

yj 10 (1998) - 164737 , i.e., a filter in which its response to 97% or 

pi more of the spatial frequency component at the grid pitch of the 

Olio stationary grid 4 is reduced to 5% or less. It can be realized, 

for example, by a wavelet transform filter employing filter 
p coefficients, such as (17, 7) - tap coefficients listed in Table 1, 

fU (13, 7) -tap coefficients listed in Table 2, and (15, 5) -tap 

O coefficients listed in Table 3. The frequency response 

15 characteristics of the low-pass filters in wavelet transform 

filters are shown in Fig. 8. 

Table 1: (17, 7) -tap 
Wavelet split low-pass filter Gl 
20 (a8,a7,a6,a5,a4,a3,a2,al,a0,al,a2,a3,a4,a5,a6,a7,a8) 
Wavelet split high-pass filter HI 
(b3 , b2 , bl, bO , bl , b2 , b3 ) 



aO = 


0 , 


. 53534743 


bO = 


1 . 


26885769 


al = 


0 . 


. 31296897 


bl = 


-0 


. 50000000 


a2 = 


= -0 


. 02883059 


b2 = 


- 0 


. 13442884 


a3 = 


= -0 


. 09095774 


b3 = 


- 0 


. 13442884 
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a4 = 0.01828153 
a5 = 0.03402453 
a6 = -0.01037186 
a7 = -0.00278855 
5 a8 = -0.00278855 

Wavelet synthesis low-pass filter Gl ' 
(-b3,b2, -bl,bO, -bl,b2, -b3) 
£i Wavelet synthesis high-pass filter HI' 

Cn (a8,-a7,a6,-a5,a4,-a3,a2,-al,a0,-al,a2,-a3,a4, 

U1 

pllO -a5, a6 , -a7 , a8) 

nj 

O Table 2 : (13, 7) -tap 

fy Wavelet split low-pass filter Gl 



(a6,a5,a4,a3,a2,al,a0,al,a2,a3,a4,a5,a6) 
15 Wavelet split high-pass filter HI 

(b3,b2,bl,b0,bl,b2 , b3) 
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aO = 


0 . 


5590000000 


bO = 


1 . 


1678751323 


al = 


0 . 


3070000000 


bl = 


- 0 


. 5408349230 


a2 = 


- 0 


. 0460000000 


b2 = 


- 0 


. 0839375665 


a3 = 


-0 


. 0901544756 


b3 = 


0 . 


0408349233 


a4 = 


0 . 


0338214681 








a5 = 


0 . 


0308330075 








a6 = 


-0 


.0150000000 









Wavelet synthesis low-pass filter Gl ' 
2 5 {-b3,b2, -bl,bO, -bl,b2, -b3) 

Wavelet synthesis high-pass filter HI' 
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(a6, -a5,a4, -a3,a2, - al,aO, -al,a2, -a3,a4, -a5,a6) 



mio 
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Table 3: (15, 5) -tap 

Wavelet split low-pass filter Gl 

(a7,a6,a5,a4,a3,a2,al,a0,al,a2,a3,a4,a5,a6,a7) 
Wavelet split high-pass filter HI 

(b2,bl,b0,bl,b2) 



aO = 


0 . 


. 53534743 


bO = 


1 , 


.26885769 


al = 


0 , 


. 31296897 


bl = 


- - 0 


. 50000000 


a2 = 


= - 0 


. 02883059 


b2 = 


- 0 


. 13442884 


a3 = 


= - 0 


. 09095774 








a4 = 


0 . 


. 01828153 








a5 = 


0 . 


. 03402453 








a6 = 


= - 0 


. 01037186 








a7 = 


- 0 


. 00278855 
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Wavelet synthesis low-pass filter Gl ' 

(b2, -bl,b0, -bl,b2) 
Wavelet synthesis high-pass filter HI' 

(-a7,a6, -a5,a4, -a3,a2, -al,aO, -al,a2, -a3,a4, 
- a5 , a6 , - a7 ) 
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With the wavelet transform filters, when a horizontal 
grid is used as the stationary grid 4, the grid component develops 
in the signal LHl . When a vertical grid is used, the grid component 
develops in the signal HLl . When a cross grid is used, the grid 
component develops in the signal HHl . In the signal LLl, a moire 
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pattern resulting from the stationary grid 4 hardly develops 
regardless of grid direction, because the spatial frequency 
component of the stationary grid 4 has sufficiently been 
suppressed . 

Among the wavelet - transform coefficient signals HHl, 
HLl, LHl, andLLl, the signal s HHl , HLl, and LHl having a possibility 
of containing the grid component are input to the direction judging 
means 35 and the grid component suppressing means 36. 

Furthermore, in the wavelet transform means 32a of the 
gilO second stage, the signal LLl is subjected to wavelet transformation 

il by the use of basic wavelet functions HO, GO, whereby signals HH2, 

B 

O HL2, LH2, and LL2 are obtained. The signal LL2 represents a 

one - sixteenth reduction image obtained by reducing the original 
image to one- fourth in both the vertical and horizontal directions . 
15 In the one - sixteenth reduction images of the original image, the 

signals HL2, LH2, and HH2 represent the images of the vertical edge , 
horizontal edge, and diagonal edge components, respectively. 
Since the grid component hardly develops in the signal LLl, as 
described above, the wavelet functions HO and GO which are employed 
20 at the second stage do not need to be set at the grid pitches of 

the stationary grid 4, unlike the wavelet functions HI and Gl 
employed at the initial stage. For example, the wavelet functions 
HO , GO may be a wavelet transform filter employing Daubechie ' s ( 9 , 
7) -tap filter coefficients listed in Table 4. The frequency 
25 response characteristic of the low-pass filter in this (9, 7)- 

tap wavelet transform filter is shown in Fig. 8. 
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Table 4: Daubechie's (9, 7) -tap 
Wavelet split low-pass filter GO 

(a4,a3,a2,al,a0,al,a2,a3,a4) 
Wavelet split high-pass filter HO 

(b3,b2,bl,b0,bl,b2,b3) 



aO 


= 0.602949180 


bO 


= 1.115087052000 


al 


= 0.266864120 


bl 


= - 0 . 591271763000 


a2 


= -0.078223267 


b2 


= -0.057543526000 


a3 


= -0.016864118 


b3 


= -0.091271763114 


a4 


= 0.026748757 






Wavelet 


synthesis low-pass filter GO' 


(-b3,b2, -bl,bO, -bl,b2 


, -b3) 




Wavelet 


synthesis high- 


pass 


filter HO' 



(a4 , -a3 , a2 , -al , aO , -al , a2 , - a3 , a4 ) 

In the same manner as the aforementioned second stage, 
a wavelet - transform coefficient signal LLk obtained for each 
frequency band is subjected to wavelet transformation by n times, 
whereby wavelet- transform coefficient signals HHl to HHn, HLl to 
HLn, LHl to LHn, and LLl to LLn are obtained . The wavelet- transform 
coefficient signals HHn, HLn, LHn, and LLn, obtained by the nth 
wavelet transformation, each represent a 1/ (1/2) ^" reduction image 
in which the number of pixels in each of the horizontal and vertical 
directions is (1/2)", compared with the original image signal S^^^. 
The greater the "n" in the wavelet - transform coefficient signals 
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HHn, HLn, LHn , and LLn, the lower the frequency band. Thus, each 
of the wavelet - transform coefficient signals HHk, HLk, LHk, and 
LLk (where k is an integer of 1 to n, representing resolution level) 
becomes a band- limited image signal carrying a frequency component 
5 which has a predetermined frequency range contained in the 

frequency range of the original image signal S^^.^. The signal HHk 
represents changes in frequency in the horizontal and vertical 
\0 directions of the original image signal S^^g, and becomes a lower 

on frequency signal if k is greater. The signal HLk represents a 

01lO change in frequency in the horizontal direction of the original 

image signal S^^^g, and becomes a lower frequency signal if k is 
greater. The signal LHk represents a change in frequency in the 
fIJ vertical direction of the original image signal S^^g, and becomes 

nu 

S a lower frequency signal if k is greater. 

15 Next, in the direction judging means 35, the power of 

each of the signals LHl, HLl, and HHl is calculated. Then, it is 
judged whether or not each power is greater than a predetermined 
threshold THO. Based on the result, the grid direction of the 
stationary grid 4 is judged . More specifically, as described above, 
20 the grid component develops in the signal LHl when a horizontal 

grid is used, develops in the signal HLl when a vertical grid is 
used, and develops in the signal HHl when a cross grid is used. 
Therefore, when only the power of the signal LH is greater than 
the threshold value THO it is judged that a horizontal grid has 
25 been used, when only the power of the signal HL is greater than 

the threshold value THO it is judged that a vertical grid has been 
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used, and when both the power of the signal HL and the power of 
the signal LH are greater than the threshold value THO it is judged 
that a cross grid has been used. This result of judgement is input 
to the stationary grid - component suppressing means 36. 

Next, based on the grid direction judged by the direction 
judging means 35, in the stationary grid- component suppressing 
means 36 a signal of the signals LHl, HLl, and HHl which contains 
a stationary grid component is further subjected to one - 
dimensional wavelet transformation recursively and repeatedly by 
a predetermined number of times in the grid length direction of 
the stationary grid 4 (grid direction) by the use of a predetermined 
band splitting filter. For instance, in the case where a vertical 
grid is used as shown in Fig. 2, the signal HLl is subjected to 
one - dimensional wavelet transformation in the vertical scanning 
direction, whereby a signal HLl (LI) , carrying a low frequency side, 
and a signal HLl (HI) , carrying a high frequency side, are obtained. 
Then, the signal HLl (LI) carrying a low frequency side is subjected 
to one - dimensional wavelet transformation in the vertical scanning 
direction, whereby a signal HLl (LI) , carrying a low frequency side, 
and a signal HLl (HI) , carrying a high frequency side, are obtained. 
Similarly, by repeatedly performing one - dimensional wavelet 
transformation on the signal HLl (Lj ) , which carries a low frequency 
side (where j is an integer of 1 to m) , in the vertical scanning 
direction by m times , the signal HLl is subj ected to one - dimensional 
wavelet transformation in the vertical scanning direction over a 
plurality of stages (m stages), whereby wavelet - trans form 



coefficient signals HLl(Hl) to HLl (Hm) , HLl (Lm) are obtained. 

Fig, 9 shows the decomposed components of the original 
image signal obtained after the original image signal has been 
subjected to wavelet transformation. In Fig. 9A, the first 
two-dimensional wavelet transformation is performed, and in Fig. 
9B, after the second two-dimensional wavelet transformation the 
signal HLl is subjected to one - dimensional wavelet transformation 
in the vertical direction once, 
gl The effects obtained by repeating one - dimensional 

Olio wavelet transformation in the grid length direction are 

illustrated in Fig. 10. 

The wavelet - transform coefficient signal HLl represents 
a low frequency component in the vertical scanning direction and, 
as shown in Fig. 10, contains both a vertical grid component and 
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py 

15 an image signal which has the vertical direction of an original 



image signal as its main component. If the signal HLl is subjected 
to one-dimensional wavelet transformation in a grid length 
direction (in this embodiment, vertical scanning direction) , the 
grid component in the case of a vertical grid is considered to be 

20 a one - dimensional frequency signal which is vertically long, and 

therefore, band split is made so that a low frequency component 
containing the vertical grid component is contained in the signal 
HLl(Lj) and the other high frequency components are contained in 
the signal HLl(Hj). If this one - dimensional wavelet 

25 transformation is repeatedly given, a low frequency signal HLl (Lm) 

in which band split has sufficiently been made represents an 
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extremely low frequency component which contains a grid component 
in the grid length direction independently of grid pitch size. If 
the "m" in the signal HLl(Lm) becomes greater, the signal will 
represent only a grid component. 

Actually, in most cases, the stationary grid is not 
perfectly placed horizontally, vertically, or in a crossing 
direction where "horizontal" and "vertical" are 1 : 1. Therefore, 
too great a number (m) of repeats will cause the grid component 
to develop in the signal HLl(Hj). Thus, it is preferable that 
one - dimensional wavelet transformation be repeated a few times. 

Next, among the signals HLl (HI) to HLl (Hm) , and HLl (Lm) 
obtained by applying one-dimensional wavelet transformation to the 
signal HLl in the vertical scanning direction over a plurality of 
stages, the signal HLl(Lm) which is a component on the lowest 
frequency side is made zero. That is, the signal HLl(Lm) 
containing a grid component is suppressed. Thereafter, the signal 
HLl(Lm), made zero, and the signals HLl (HI) to HLl (Hm) , are 
subjected to inverse one - dimens ional wavelet transformation, 
whereby a signal HLl' is obtained. As previously described, at 
least the signal HLl (Lm) containing the spatial frequency 
component of the stationary grid 4 has been suppressed. Therefore, 
the signal HLl', in which the spatial frequency component with a 
predetermined frequency range containing the grid component has 
been reduced, is obtained. 

Next, in the processing means 33 a desired process (e.g. , 
an enhancing process) is applied as occasion demands . Thereafter, 





in the inverse wavelet transform section 34, the signals LLn, HLk, 
and LHk are sequentially subjected to inverse wavelet 
transformation from level n to level 1. 

Fig. 11 shows the construction of the inverse wavelet 
transform section 34. As shown in the figure, in the inverse 
wavelet transform means 34a the lowest frequency band signals HHn, 
HLn, LHn, and LLn are subjected to inverse wavelet transformation, 
whereby an signal LLn-1 is obtained. 



performed by each wavelet transform means 34a. As shown in the 
figure, signals LLn (LLk) and LHn(LHk) are subjected to a process 
for leaving space for 1 pixel between adjacent pixels, in the 



12. The signals are also subjected to a filtering process in the 
vertical scanning direction by the inverse wavelet transform 
functions GO', HO' corresponding to the functions GO ' , HO' employed 
in performing wavelet transformation , and are added . Furthermore , 
the signal obtained by addition (referred to as a first added 
signal) is subjected to the aforementioned space leaving process 
in the horizontal scanning direction, and is subjected to the 
filtering process in the horizontal direction by the function GO ' . 
In this manner, a first signal is obtained. On the other hand, 
signals HLn (HLk) and HHn(HHk) are subjected to the process for 
leaving space for 1 pixel between adjacent pixels, in the vertical 
direction. The signals are also subjected to the filtering process 
in the vertical scanning direction by the inverse wavelet transform 



Fig. 12 shows how the inverse wavelet transformation i 



s 



vertical direction. 



This process is represented by " t2" in Fig 



functions GO', HO', and are added. Furthermore, the signal 
obtained by addition (referred to as a second added signal) is 
subjected to the process for leaving space for 1 pixel between 
adjacent pixels, in the horizontal scanning direction, and is also 
subjected to the filtering process in the horizontal direction by 
the function HO'. In this way, a second signal is obtained. The 
first signal and the second signal are added, whereby a signal 
LLn-l(LLk-l) is obtained. 

Next, in the inverse wavelet transform means 34a, signals 
HHn-1, HLn-1, and LLn-1 are subjected to inverse wavelet 
transformation in the same manner as the aforementioned, whereby 
a processed signal LLn-2 is obtained. In the same way as the 
aforementioned, a signal LLl is obtained by repeating inverse 
wavelet transformation down to resolution level 1. 

At the resolution level 0 representing the original image, 
the signals LLl , LHl , HHl , and HLl ' are subj ected to inverse wavelet 
transformation by the inverse wavelet transform functions Gl ' , HI ' 
corresponding to the functions Gl, HI employed in performing 
wavelet transformation at the initial stage, whereby an image is 
reconstructed. Note that the signal HLl' has been reduced in grid 
component. In the reconstructed image, therefore, the grid 
component of the stationary grid 4 has been reduced. 

Fig. 13 is used to explain advantages of the present 
invention. Fig. 13A illustrates a frequency response 
characteristic obtained in the case of a method described in 
Japanese Unexamined Patent Publication No. 1 0 ( 1 9 9 8 ) - 1 64 7 3 7 
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(hereinafter referred to as a conventional method). Fig. 13B 
illustrates a frequency response characteristic obtained in the 
case of the present invention, and Fig. 13C illustrates the 
frequency response characteristics in Fourier space obtained by 
the conventional method and the method according to the present 
invention. Note that Fig. 13C illustrates the case where a 
vertical stationary grid is employed in the Fourier space with the 
horizontal scanning direction as a v-axis and the vertical scanning 
direction as a u-axis. 

In the case of the conventional method, a filter is 
employed in which not only its response at a frequency near a spatial 
frequency corresponding to the grid pitch of a stationary grid, 
but also its response at a high frequency component greater than 
that, is made zero. Therefore, as shown in Fig. 13A, the high 
frequency component greater than the spatial frequency 
corresponding to the grid pitch is also suppressed. As a result, 
not only the stripe pattern of the stationary grid 4, but also the 
high frequency component originally contained in the image, is 
removed and reduced. Because of this, the image will be reduced 
in sharpness. In addition, because the filtering process in the 
conventional method is not a filtering process taking direction 
into consideration, all (oblique line portion) other than a low 
frequency domain (central blank portion) on the Fourier space is 
suppressed as shown Fig. 13C. That is, not only a vertical pattern 
due to a vertical grid (which is desired to be really suppressed) 
but also a high frequency component, which is not desired to be 



suppressed, such as a horizontal pattern, a diagonal pattern, etc. , 
which are contained in the original image, is suppressed. 

On the other hand, in the case of the present invention, 
even if a filter exhibiting the same characteristic as that 
5 described in the aforementioned Japanese Unexamined Patent 

Publication No. 10 (1998) -164737 is used as the filter employed in 
the wavelet transformation at the initial stage, only a 
yi predetermined range near the spatial frequency component of the 

m stationary grid 4 can be suppressed (or cut) and the remaining high 

U1 

n^lO frequency components sustained, as shown in Fig. 13B. Therefore, 

ru 

L-L an image, in which a stripe pattern due to the stationary grid is 

p inconspicuous and sharpness is high, can be obtained . In addition, 

f|j the signals HLl, LHl, and HHl can be switched according to the grid 

£2 direction so that only a signal of the signals HLl, LHl, and HHl 

15 which contains the grid component is suppressed. Therefore, as 

shown in Fig. 13C, only a predetermined range on the high frequency 
side can be suppressed. The predetermined range on the high 
frequency side is in a domain (with a slight width in a direction 
perpendicular to grid direction) near the grid- direction axis (in 
20 the case of a vertical grid, v-axis) of the Fourier space, and 

contains the spatial frequency component of the stationary grid 
4. For example, in the case of a vertical grid, only a vertical 
pattern due to the grid is suppressed and there is no possibility 
that a high frequency component, such as a horizontal pattern, a 
25 diagonal pattern, etc., which is contained in the original image, 

will be suppressed. 
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In addition, the direction of a grid can be judged, based 
on the wavelet - transform coefficient signals HHl, HLl, and LHl 
having a possibility of containing the grid component, among the 
wavelet - transform coefficient signals HHl, HLl, LHl, and LLl 
obtained by two-dimensional wavelet transformation at the initial 
stage. Therefore, when performing the process of suppressing a 
grid component, there is no necessity for previously knowing what 
kind of grid is used in photographing. 

Furthermore, if signals are restored to resolution level 
1 without being restored to resolution level 0, a stripe pattern 
resulting from the stationary grid 4 will barely develop, because, 
as previously described, in the signal LLl the spatial frequency 
component of the stationary grid 4 is sufficiently suppressed. 
Moreover, since the signal LLk after level 1, obtained by giving 
wavelet transformation to the signal LLl in which the grid component 
has been suppressed, does not contain a moire component, there is 
no possibility that in all reduced - scale images, a moire pattern 
resulting from the stationary grid 4 will occur. Even if the image 
is arbitrarily enlarged or reduced, no moire pattern will occur. 
Therefore, even if an image is subjected to a desired process such 
as an enhancing process and the processed image is restored by 
inverse mul tiresolution transformation, the restored image will 
contain no artifact due to the moire pattern even if it has any 
resolution level. As a result, a high-quality image which is easy 
to view can be provided. 

Thus, since the grid component is suppressed by wavelet 




transformation (mul tiresolution transformation) at the initial 
stage, a moire pattern will no longer develop in the images having 
the subsequent resolution levels. In addition, because wavelet 
transformation is employed, this embodiment is suitable for 
5 obtaining an enlarged- or reduced - scale image and is also 

convenient in transferring image signals through networks. 

While the present invention has been described with 
J4 reference to a preferred embodiment thereof, the invention is not 

^ to be limited to the details given herein, but may be modified within 

Si 10 the scope of the invention hereinafter claimed. 

m 

For example, while it has been described that the 
stationary grid 4 is a vertical grid, if it is other than the 
vertical grid the signal HLl to be handled may be the signal LHl 
or HHl, depending on the grid direction. For example, in the case 

I It 

15 where the stationary grid 4 is a horizontal grid, the signal LHl 

is further subjected to one - dimensional wavelet transformation in 
the horizontal scanning direction over a plurality of stages. 
After the signal component on the low frequency side is made zero, 
inverse one - dimensional wavelet transformation is applied to 

20 reconstruct the image. In the case where the stationary grid 4 

is a cross grid, the signals HL and LH are further subjected to 
one - dimensional wavelet transformation in the vertical and 
horizontal scanning directions over a plurality of stages, 
respectively. After the signal component on the low frequency side 

25 is made zero, inverse one - dimensional wavelet transformation is 

applied to reconstruct the image. 
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Instead of applying one - dimensional wavelet 
transformation to the signal HL and/or signal LH, a high-pass filter 
for reducing a grid component may be employed. 

If a stationary grid to be used is subjected to the 
aforementioned processing in each grid direction independently of 
the direction of a grid which is actually used, for example, and 
if in addition to the processing shown in Fig. 9, the signal LHl 
is repeatedly subjected to one - dimensional wavelet transformation 
in the horizontal scanning direction, the suppressing effect can 
be obtained not only in a predetermined range on the high frequency 
side near the v-axis shown in Fig. 13C, which contains the spatial 
frequency component of the stationary grid 4, but also in a 
predetermined range on the high frequency side near the u-axis, 
which contains the spatial frequency component of the stationary 
grid 4. In many cases, a grid of the same direction as either of 
horizontal and vertical scanning directions is usually employed. 
In those cases, the grid component will be contained in either the 
signal HL (when a vertical grid is used) or signal LHl (when a 
horizontal grid is used) , and consequently, the effect of the 
aforementioned suppressing process can be obtained without 
considering grid direction. Also, in those cases, since the 
suppressing process is not applied in infinite directions but is 
limited only to each grid direction in which the suppressing process 
is to be performed, the influence of the suppressing process on 
an image can be reduced, compared with prior art . Note that a cross 
grid can be handled by repeatedly applying one - dimensional wavelet 



transformation to the signal HHL in both the horizontal and vertical 
scanning directions . 

In this embodiment, it has been described that the image 
signal S is subjected to two-dimensional wavelet transformation. 
5 However, in the case where grid direction is known, the original 

image signal may be subjected to one - dimensional wavelet 
transformation in that direction (in the vertical scanning 
direction when the grid is a vertical grid) , and the wavelet- 



y3 

x| transform coefficient signal containing a grid component (signal 

Un 10 Llof Fig. 14B) sub j ec ted to the process of reducing a grid component , 

ny as shown in Fig. 14. Note that the wavelet function which is 



}""■ 

s employed in performing one - dimens ional wavelet transformation may 

Q 

be an arbitrary one. For instance, the aforementioned functions 

m 

nJ HO, GO can be employed. 

15 If the method shown in Fig. 14 is employed, an extremely 

low frequency component representing grid- length direction will 
develop in a low frequency wavelet - transform coefficient signal 
(signal LI) of the wavelet - transform coefficient signals. After 
the wavelet- transform coefficient signal LI has repeatedly been 

20 subjected to the process of reducing a grid component (e.g., 

one - dimensional wavelet transformation in the same direction) (see 
Fig. 14C) , the lowest frequency wavelet - transform coefficient 
signal is made zero, or only signal components representing the 
original image other than the grid component are extracted with 

25 a high-pass filter. In this way, an extremely low frequency 

component representing the grid component can be suppressed. 
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Thus, when an image is restored by applying one- 
dimensional wavelet transformation to the signal in which the grid 
component has been reduced, an image with a reduced grid component 
can be restored regardless of the resolution level at which the 
5 image is restored . 

In the embodiment shown in Fig. 14, the above-mentioned 
processing can also be applied in each grid direction of a 
stationary grid to be used, independently of the grid direction 
that is actually used. For example, in addition to the processing 
10 in the vertical scanning direction shown in Fig. 14, the processing 

is performed in the horizontal scanning direction. Between two 
images finally obtained, the image with fewer grid components can 
be used. 

It has been described that when the image signal S 
15 represented in a real space domain is transformed into a plurality 

of image signal which can be handled in a frequency domain, the 
mul tiresolution decomposing process employing the wavelet 
transform process is used. However, any method of transformation 
can be employed as long as it can apply a process of reducing a 
20 component which has a desired frequency range containing a 

stationary grid component contained in an image signal. For 
example, a mul tiresolution decomposing process employing a 
Laplacian pyramid expansion can also be employed. Also, an image 
signal expressed in the real space domain may be transformed by 
25 Fourier transformation into image signals expressed in a frequency 

domain (frequency spectra). 
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In addition, all of the contents of Japanese Patent 
Application Nos. 2000-011174 and 2000-395577 are incorporated into 
this specification by reference. 

5 
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